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Abstract 

A  model  for  the  underpotential  deposition  of  metals  that  occurs  in  stages  is  introduced. 
In  this  model  the  deposition  takes  place  as  a  sequence  of  first  order  phase  transitions  of  the 
adsorbate.  In  this  application  we  study  the  underpotential  deposition  of  Cu  on  a  Au(lll) 
surface  in  the  presence  of  sulfate  ions.  The  voltammogram  of  the  deposition  shows  two  sharp 
spikes  which  are  reproduced  by  our  model. 
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1  INTRODUCTION 


\  oltanmiograms1  of  underpotential  deposited  films  on  perfect  single  crystal  surfaces 
contain  sharp  spikes.2-6  Such  spikes  seem  to  occur  also  in  polymer  film  formation.'  *  In  earlier 
work,  we  discussed  the  possibility  of  explaining  these  spikes  as  the  result  of  first  order  phase 
transitions  occuring  in  the  surface.1  10  There  are  a  number  of  conditions  that  have  to  be  met 
to  obtain  sharp  spikes  in  the  voltammogram.  These  include  chemical  equilibrium,  the  degree 
of  perfection  of  the  substrate  (a  single  crystal  in  most  cases),  and  the  scanning  rate  of  the 
voltammogram.  Ideally  the  width  of  the  voltammogram  is  due  to  diffusion  processes.1  so 
that  slow  scanning  rates  would  be  best  to  observe  these  spikes.  The  area  under  the  spike  is 
proportional  to  the  change  in  coverage  of  the  surface,  which  is  not  necessarily  equal  to  the 
ideal  coverage.  As  was  discussed  in  earlier  work.9  "  second  order  phase  transitions,  such  as 
the  one  occuring  in  the  hard  hexagon  lattice  gas  model,  are  also  possible.  However,  these 
transitions  would  be  very  hard  to  see  in  a  voltammogram.  because  of  the  vanishingly  small 
area  of  the  spike  associated  with  a  second  order  transition. 

We  have  recently  studied  a  statistical  mechanical  model  for  adsorption  on  electrode 
surfaces  that  includes  the  lateral  interactions  of  the  adsorbates.9,10  In  this  model  we  use 
exact  results  for  two  dimensional  lattices  that  tell  us  when  a  first  or  a  second  order  phase 
transition  will  occur.  We  also  show  that  a  first  order  phase  transition  corresponds  to  a  sharp 
spike  in  the  voltammogram.10 

In  some  systems,  more  than  one  spike  has  been  observed  to  occur  in  the  voltammogram 
during  the  deposition  process.  In  particular,  a  voltammogram  with  two  sharp  spikes  has 
been  observed  during  the  underpotential  deposition  of  copper  on  the  (111)  surface  of  gold 
in  the  presence  of  sulfate  ions.2  This  system  has  been  recently  studied  by  a  number  of 
in-situ  techniques.1216  While  we  do  not  have  a  complete  and  clear  picture  yet.  the  recent 
scanning  tunnel  microscope  and  atomic  force  microscope  experiments  seem  to  indicate  that 
there  is  a  \/3  x  \/Z  structure  that  exists  between  the  two  underpotential  deposition  peaks 
at  0.07  and  0.22  volts.  There  seems  to  be  evidence  also  that  there  is  sultate  adsoiptioii.1 ' 
which  undoubtedly  plays  an  important  role  in  the  underpotential  deposition.  Unfortunately, 
careful  measurements  of  the  fraction  of  occupied  surface  are  not  at  hand.  A  preliminary 
measurement  of  the  relative  area  of  the  two  spikes  of  the  voltammogram  seems  to  show 
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that  this  ratio  is  1  to  2. 18  We  have  constructed  a  physically  reasonable  mode)  which  is 
consistent  with  the  near  edge  EXAFS.11  and  which  exhibits  two  first  order  phase  transitions 
corresponding  to  a  voltammogram  with  two  spikes  of  this  relative  area.  For  a  reasonable 
choice  of  parameters,  it  can  be  made  to  agree  with  the  experimentally  obtained  one.  Treating 
incommensurate  structures  such  as  the  ones  recently  postulated12  n  is  certainly  possible  in 
our  theory,  but  we  could  not  use  the  simple  and  exact  results  of  two  dimensional  lattices.  Our 
model  is  simple  and  versatile,  and  can  be  adapted  to  other  circumstances,  such  as  different 
lattices  or  adsorption  sequences. 


2  THE  MODEL 

In  our  model  for  the  underpotential  deposition  of  Cu  on  Au(lll)  in  the  presence  of 
high  concentrations  of  sulfuric  acid  we  assume  that  a  well  defined  sequence  of  events  takes 
place: 

•  First. the  sulfate  ions  that  are  adsorbed  form  a  \/3  x  \/3  lattice  on  the  gold  surface. 

•  Copper  ions  are  then  adsorbed  on  the  free  adsorption  sites,  which  form  a  honeycomb 
lattice. 

•  Finally  the  adsorbed  sulfate  ions  are  replaced  by  copper  ions  . 

We  assume  that  the  sulfate  binds  to  the  gold  (111)  surface  in  such  a  manner  that  the 
sulfur  is  directly  on  top  of  the  adsorption  site  for  the  copper,  three  of  the  sulfate  oxygens 
being  above  and  directly  associated  with  the  three  gold  atoms  of  the  surface,  which  form  a 
triangle  about  the  adsorption  site.  Packing  considerations  indicate  that  two  sulfate  groups 
cannot  be  adsorbed  onto  neighboring  adsorption  sites.  The  sulfate  ions  will  thus  form  a 
\/3  x  \/3  film  by  occupying  one  of  the  three  triangular  sublattices  Ar  of  the  full  triangular 
lattice  of  adsorption  sites.  One  can  show  from  statistical  mechanics  that  an  order-disorder 
transition  will  accompany  this  adsorption  (or  desorption)  process.  However,  since  then'  is 
no  change  in  the  adsorption  coverage  at  the  transition  (  rather  a  change  in  the  adsorption 
rate  ).  there  should  be  no  observable  spike  in  the  voltammogram.  This  type  *  /I 

sulfate  film  is  illustrated  in  Fig.  1.  The  gold  atoms  are  represented  by  large  white  disks, 
the  adsorption  sites  by  small  black  disks.  Each  adsorbed  sulfate  ion  is  represented  by  a  set 
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of  three  lines  centered  at  the  adsorption  site  and  ending  at  the  three  associated  gold  atoms. 
Each  line  represents  a  S-0  bond.  The  vacant  adsorption  sites  form  a  honeycomb  lattice  \H. 

As  the  potential  bias  becomes  less  positive,  copper  atoms  will  be  adsorbed  onto  the 
surface  and  discharged.  To  model  this  deposition  process,  we  construct  a  statistical  mechan¬ 
ical  model  for  the  interaction  of  a  solution  of  copper  ions  with  a  gold  (111)  surface  which 
contains  the  \/3  x  \/3  sulfate  film. 

As  was  done  before  ,  the  adsorption  sites  in  the  model  are  considered  to  be  stickv 
points  (to  which  atoms  are  either  chemically  bound  or  not)  with  a  stickiness  parameter  A  7- 
for  the  sites  on  A j.  the  triangular  sublattice  of  adsorption  sites  which  are  associated  with 
the  sulfate  groups,  and  with  a  stickiness  parameter  A//  for  sites  on  the  vacant  honeycomb 
sublattice  A#. 

This  sticky  site  model  can  be  shown  to  be  equivalent  to  a  two  dimensional  lattice 
gas 9.19.20  Qf  COpper  atoms  which  have  a  fugacity 

*r  =  *TPi(0.\ *)  (1) 

if  adsorbed  on  Aj.  and  a  fugacity 


=  W”(0.  $)  (2) 

if  adsorbed  on  A//.  Here  p°(0.  'f)  is  the  number  density  of  copper  ions  at  contact  with  the 
surface  under  a  potential  bias  Here 

*  =  3([v(0)  -  vpzc).  (3) 

where  e  is  the  elementary  charge.  3  =  1/kT,  c  (0)  is  the  potential  at  the  electrode  surface, 
and  u'p:c  is  the  potential  of  zero  charge.  At  23°C 

=  38.9221*  (4) 

where  V  is  the  potential  in  volts.  We  approximate  p°( 0.  $)  as10 

/>?((>,*)  =  f-J|V(0.0).  (5) 
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where  £1  =  2  is  the  electrovalence  of  the  copper  ions  in  solution.  When  two  copper  ions  are 
adsorbed  on  sites  i  and  j.  separated  by  a  distance  ru.  they  interact  in  the  equivalent  lattice 
gas  with  an  interaction  energy  which  is  equal  to  the  potential  of  mean  force 

w(rtj)  =  -kT\n\s,°i(rij)].  (6) 

where  is  the  pair  correlation  function  of  two  copper  atoms  in  the  adsorption  plane.  We 

should  remark  that  this  pair  correlation  function  is  quite  difterent  from  the  pair  correlation 
function  of  the  copper  ions  in  the  bulk,  since  they  already  form  part  of  the  metal  electrode. 

We  shall  calculate  the  properties  of  the  two  dimensional  lattice  gas  using  the  following 
simplifying  assumptions:  if  two  copper  atoms  are  adsorbed  on  neighboring  sites  of  A//,  their 
interaction  energy  is 

wH  =  -k'Tln[g°2(d)l  (7) 

and  if  they  are  adsorbed  on  neighboring  sites  of  Aj.  then  they  interact  with  an  energy 

tty  =  -ATln[</°(^3d)].  (8) 

Otherwise  we  shall  assume  that  there  is  no  interaction  between  adsorbates.  This  approxi¬ 
mation  decouples  the  adsorption  process  on  the  two  sublattices,  and  consequently  simplifies 
the  calculation.  (  We  will  include  the  interactions  between  atoms  adsorbed  on  \j  and  atoms 
adsorbed  on  A h  in  a  forthcoming  publication.)  . 

We  shall  now  use  some  ideas  introduced  in  Ref.(  10)  to  calculate  the  adsorption  isotherms 
for  our  model.  We  use  the  first  8  terms  in  the  low  zj  series  approximation  to  the  free 
energy21  22  to  construct  a  Pade  approximant  to  9j,  the  fraction  of  sites  of  A 7-  which  are 
occupied  by  copper  atoms.  This  Pade  approximant.  9j is  best  at  low  zj  and  is  constructed 
to  be  the  natural  extension  of  the  Langmuir  isotherm  which  includes  interactions  between 
atoms  adsorbed  on  neighboring  sites.  Using  symmetry  considerations  we  then  costruct  a 
Pade  approximant  9t,k  which  is  best  at  high  zj.  The  explicit  form  of  these  Pade  approxi¬ 
mant  s  is  given  in  Ref.  (10). 

When  the  lateral  interactions  are  turned  off.  </2(v/3</)  =  1  and  both  9T  h  and  0t.i  are 
equal  to  the  exact  isotherm.  9j,  which  is.  in  this  case,  identical  with  Langmuir's  isotherm. 


The  isotherm  Bj  is  a  smooth  function  of  zj  for  all  values  of  g°(\/3d)  which  are  less  than 
the  critical  value  3. 9  For  values  exceeding  3.  there  is  a  first  order  phase  transition  with  a 
corresponding  jump  in  the  fraction  of  occupied  sites  Bj  . 

In  order  to  construct  a  single  approximation  to  Bj  which  is  valid  for  all  values  of  zj. 
we  use  a  switching  function 

^(</r--sr)  =  ( 1/2) {1  +  Erf[(yT  -  1).S7]}.  (9) 

where  Erf  is  the  error  function,  and 

VT  =  zT-lg°2(V3d)~:i  (10) 

is  a  variable  which  becomes  unity  at  the  first  order  phase  transition.9  The  parameter  sj  is  a 
measure  of  the  sharpness  of  the  switching  function  and  is  introduced  in  order  to  account  for 
a  variety  of  line  broadening  factors,  especially  diffusion  and  the  degree  of  crystallinity  of  the 
surface,  which  are  not  included  explicitly  in  our  present  treatment  of  the  model.  We  then 
construct  the  single  continuous  approximation 


=  g(yr- $t)Bt.(  +  [1  —  yiyr-^T^T.h-  (11) 

In  a  similar  fashion  we  use  8  terms  in  the  low  zw  approximation  to  the  free  energy  of 
a  lattice  gas  on  the  honeycomb  lattice  A#  to  construct  a  Pade  approximant  BH ,  which  is 
best  at  low  and  a  Pade  approximant  Bfj^  which  is  best  at  high  zg.  The  fraction  of  sites 
of  \n  occupied  by  copper  atoms  is  then  approximated  as 


&H  =  y(yH**H)0H.t  +  [1  -  'l{yH-*H)]&H.h-  (12) 

where 

yH  =  ZH-'gt(<irV2  (13) 

is  a  variable  which  is  unity  at  the  phase  transition  of  A#. 

The  fractional  coverage  of  the  full  lattice  A 7  U  A h  is  then  given  as 

B  =  (2/3)0//  +  (1/3)07 .  (14) 
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The  change  in  6  at  a  phase  transition  is  roughly  equal  to  the  measured  area  under  the 
corresponding  spike  of  the  voltammogram. 


As  was  discussed  in  our  earlier  work.10  a  phase  transition  will  appear  in  the  voltam¬ 
mogram  as  a  sharp  peak  in  the  intensity  1 .  If  the  scanning  rate  is  constant  and  we  neglect 
diffusion  and  double  layer  effects.1  the  intensity  of  the  current  in  the  voltammogram  can  be 
obtained  by  differentiation  as 


'  34  ,H 


In  Fig.  2  we  show  a  voltammogram  which  is  obtained  using 


(lo) 


Ar/>?(0.0)  =  9.1. 
,W?(0.0)  =  5.2  x  10\ 
g°2((I)  =  14. 
v/3</)  =  3.1. 


and 

■SH  =  -• 

The  intensity  is  given  in  units  of  the  sweep  rate  \<I\'/Jt\.  This  voltammogram  is  very 
similar  to  the  one  obtained  experimentally  for  the  underpotential  deposition  of  copper  on  a 
(111)  surface  of  gold  in  the  presence  of  sulfuric  acid.2 
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FIGURE  CAPTIONS 


Figure  1:  The  \/3  x  \/3  sulfate  phase  on  the  (111)  face  of  gold.  Gold  atoms  are 
represented  by  large  white  disks,  the  adsorption  sites  for  sulfate  and  copper  are  depicted  as 
small  black  disks,  and  the  adsorbed  sulfate  groups  are  depicted  as  sets  of  three  lines  emerging 
from  the  adsorption  sites  to  the  neighboring  gold  atoms. 

Figure  2:  A  voltarnmogram  with  two  sharp  peaks  corresponding  to  two  first  order 
phase  transitions  in  the  model. 
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